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Abstract 

A simple, self-contained proof is presented for the concavity of the map 
{A,B) TtApK^B^^pK. The author makes no claim to originality; this 
note gives Lieb's original argument in its simplest, rather than its most gen- 
eral, form. A sketch of the chain of implications from this result to concavity 
of ^ 1-^ gA'+iog A jg presented. An independent elementary proof is 
given for the joint convexity of the map {A, B,X) Tr A+ui -^ B+ui ^'^ 

which plays a key role in entropy inequalities. 



1 Introduction 



Properties of quantum entropy, particularly the inequality known as strong subad- 
ditivity (SSA), play an important role in quantum information theory. The original 
proof of SSA is based on a concavity result of Lieb, and several approaches to SSA 
use either this result or one of his related convex trace functions ^ . These results 
have acquired an undeserved reputation as difficult to prove; indeed, the influential 
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02-1-0065, and by the National Science Foundation under Grant DMS-0314228. 
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book by Nielsen and Chuang |T2] states on p. 645 that "no transparent proof of 
SSA is known." This note is intended to remedy this situation. 

During recent lectures, I presented Lieb's original proof of the joint concavity of 
the map {A, B) Tr A^K'^B^'^K and was reminded just how simple and elegant it 
really is. Unfortunately, some of this simplicity was lost when published in a form 
that lent itself to generalizations, such as the concavity of {A,B) Tr A^K'^B'^K 
with p + q < 1. In view of the renewed interest in this result, and frequent reference 
to the elementary, but long, proof^ given in an Appendix to Nielsen and Chuang 
|12j . it seemed worth making a simple argument available to a larger audience. 

Lieb's results in j7] include the following three theorems which we consider only 
for finite dimensional matrices. 

Theorem 1 Let K be a fixed n x n matrix. For < p < I, the map {A,B) i— >■ 
F{A, B) = Tr A^ K"^ B^~'P K is well-defined on the cone of pairs of positive, semi- 
definite n X n matrices A, B. Moreover, for each p G (0, 1), the function F{A, B) is 
a concave map from this cone to [0, oo). 

Theorem 2 The map {A, B, K) ^ F{A, B, K) = Tr K^^^K^^du is well- 
defined when A, B are positive, semi-definite n x n matrices, ker(A) C ker(i^r) and 
ker(y4) C keT{K''). Moreover, F{A, B, K) is jointly convex in A, B, K . 

Theorem 3 Let K he a fixed, self-adjoint n x n matrix. The map A i— > F{A) = 
2ygA+iog(A) well-defined on the cone of positive definite matrices. Moreover, F{A) 
is concave on this cone. 

Lieb's proof of Theorem ^ uses the fact that the modulus of a function which 
is analytic and uniformly bounded on a strip is bounded by its supremum on the 
boundary. In order to make this note self-contained and accessible to readers with 
varied background, we explain this result and sketch a proof in Appendix El The 
three functions in the theorems above satisfy a homogeneity condition of the form 
F{XA) = XF{A) or F{\A,\B) = \F{A,B), etc. which has useful consequences 
summarized in Appendix iBl 

In an earlier review the author emphasized the role of the related concavity 
of A I— > Xr e^"*"^"^"^, and presented Epstein's proof |SI of this result. In fact, both 
Lieb's and Epstein's methods can be used to prove a much larger class of inequalities 
whose equivalence was established by Lieb in |7j. Epstein gives a more direct route to 
concavity of A ^jj, gX+iogA. ^j^Qg Lieb's is simpler for {A,B) i— > Tr A^ K"^ B^~p K . 
Epstein's approach requires some deep results from complex analysis, while Lieb's 
^ Based on an argument of Uhlmann JS] as presented by Simon 17, and Wehrl |20|. 
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argument requires only the maximum modulus principle. However, in the general- 
ization to TtApR^B'^K, Lieb also uses a result about concave operator functions. 
The theory of monotone and convex operator functions is closely connected to Ep- 
stein's approach and has other applications in quantum information theory. (For 
some examples and references, see 0.) 

In ; the author showed how to use the concavity of Tr e^^^°^^ to obtain simple 
proofs of the strong subadditivity of quantum entropy and the joint convexity of 
relative entropy. The advantage to the presentation in is the ease with which 
equality conditions are obtained. However, the proof of concavity of {A, B) h-> 
Tt K"^ B^~'P K given here also yields a short route to the same group of entropy 
inequalites. By observing that 



one can see that Theorem ^ also yields a simple proof (first noted in ) of the 
joint convexity of relative entropy. 

In section 3, Lieb's arguments showing a series of implications from Theorem^ 
to Theorems 121 and El are sketched. The goal is only to give the reader a feel for the 
ideas connecting seemingly disparate results. For more details, the presentation in 
Chapter 3 of ^3] is recommended. 

In Section ^ a simple new proof of Theorem |21 is presented, together with some 
remarks about its connection to entropy inequalities. 

2 Lieb's proof of Theorem [T] 

Theorem^ will first be proved for the special case A = B. Concavity then means 
that for each pair Ai, A2 > with Ai -|- A2 = 1, 



where C = XiAi + X2A2. Observe that (</), C(/)) = ^ (0, Aicj)) = {(/), ^20) = 
so that (|21) holds trivially for (p G ker(C). Hence, it suffices to prove the inequal- 
ity on ker(C)-'" so that we can assume without loss of generality that C is in- 
vertible. In finite dimensions this also implies that is bounded. Now define 



M = C^^-p)/2kCp/^ and 

fkip) = Tr AlC-P^^M^C-^^-P^/^Al-^C-^^-P^/^MC-P/^ k = 1,2. (3) 



lim i (Tr A^-^B^ -TiA) = -Tr y4(log A - log B) 



(1) 



XiTtA^K^AI-^K + \2i:tAIk^a\-''k < TtC^K^C^-pK 



(2) 



Then ^ is equivalent to 



fip) = Ai/i(p) + A2/2(p) < Tr MtM. 



(4) 
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Observe that the functions above can be analytically continued to the strip < 
Re (z) < 1. 

The next step is to show that each fk{z) is bounded on this strip. To do this, 
write z = X + iy with x,y real, and observe that for A > 0, \\A^y\\ = 1 and \\A^\\ = 
sup {i/j, A^ifj) = (sup||^||^]^('0, A?/)))'^ = \\A\\^ is the (largest eigenvalue of Ay. 
Then by repeated application of the inequality |TrXy| < Tr < ||X||ooTr|y| 
(where ||X||oo = ||-^|| is the usual operator sup norm as above), one finds that 

\fk{z))\ < \\Ak\\ \\C-'\\TiM^M<\\Ak\\ \\C-'\\ \\C\\TtK^K. (5) 

for < X < 1. Thus the functions fk{z) are uniformly bounded on the strip 
< Re{z) < 1. 

By the maximum modulus principle f Appendix . \fk{z)\ is bounded by its 
supremum on the boundary of this strip, i.e., for z = + iy or z = 1 + iy. Now, 

fk{0 + ty) (6) 

has the form TrX^F which is bounded above by (Tr X'^'X Tr y^y) Since oper- 
ators of the form yl** are unitary for t real and A positive, one finds 

|/fc(0 + iy)\< Ti M^C'y/^C-^'^AkC~^/^C-'y/^M k = l,2. (7) 

Thus 

\fiO + iy)\ < X,\f,(0 + ty)\ + X,\f,(0 + ty)\ (8) 

< Tr M^C*^/2^-^/2 (^^^^^ ^ ^^^^^ ^-l/2^~iy/2^ 

= Ti M^c'y/^c-'y/^M = ttm^m (9) 

since C was defined as Ai^i + X2A2- One can similarly show that |/(1 + iy)\ < 
Tr M^M which implies (j3)). This establishes the concavity of A Tr A^ K''' A^~^ K . 

The general case then follows from the observation 

..™-..^..(^ °)'-(° °). ,0, 

Extension to infinite dimensions: The restriction to finite dimensional matrices 
was used only to ensure that is bounded on the orthogonal complement of 
ker(C) so that (0) gives a uniform bound for \fk{z)\ on the strip < Re (2:) < 1. It 
is worth emphasizing that in this part of the proof it is enough to show that |/fc(2;)| 
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satisfies some upper bound, which can be rather crude as long as it holds uniformly 
for all z in the infinite strip. Only for the subsequent estimate on the boundary do 
we need a precise bound of the form Q, which can be generalized to operators on 
infinite dimensional spaces. 

Therefore, the theorem can be extended to infinite dimensions in several ways. 
First, observe that C = XiAi + X2A2 implies that for < < 1 

Ak < A-^C ^ Al< A-^C" C-'^/^AlC-'^/^ < A""/ (11) 

where the first implication uses the operator monotonicity of the map A ^ A'^. Then 
under the additional hypothesis that K is Hilbert-Schmidt (i.e., Ti K'^K < 00), one 
can replace (0) by 

\fk{z)\ < X^^TtM^M < Afc^ TtK^K. (12) 

Lieb uses the even weaker assumption that M = C^'^KC^^'^ is Hilbert-Schmidt, 
and also proves concavity for the map {A, B) 1— Ti A^K'^B'^K for < p+q < 1. For 
this, he uses the operator concavity of the map A A^ for < t < 1 to conclude 
that C-'/^{XiA{ + A24)C-*/2 < /. 



3 Connecting the concavity theorems 

3.1 Non- commutative multiplication and differentiation 

The key to connecting the two concavity results mentioned above was Lieb's real- 
ization that, for any fixed positive semi-definite matrix A, the following two linear 
maps on n X n matrices are inverses of each other, i.e., 

X = nA{K) = [ APRA^-Pdp (13) 
/"°° 1 1 

This is far from obvious, but can be verified by expanding in a basis of eigenvectors of 
A. Moreover, Qa can be regarded as a non-commutative version of multiplication by 
A and as a non-commutative version of multiplication by A~^. Both operators 
are positive semi-definite with respect to the inner product TiA^B = {A,B), i.e., 
TvK^nAiK) > and TiX^n^^X) > 0. 
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The operator arises when one uses the integral representation 

1 1 



log P - log Q 



Q + uI P + uI 
Q + uI^^ P + ul 



du 
1 



du 



(15) 
(16) 



to compute derivatives which arise in studying entropy. In particular, when /(x) = 
\og{A + xK) with K = self-adjoint, it follows from (Htj) that /'(O) = VL-j^^{K) 
and /"(O) = -2Ta{K) where 



:K- 



-du. 



A + uI A + uI A + uI 
This leads |13j to the useful (and norm convergent) expansion 

\og{A + xK) = log(A) + xn^\K) - x^Ta{K) + ... 



(17) 



3.2 From Theorem [T] to Theorem [2] 

To show that Theorem ^ implies Theorem |2l first observe that the joint convexity 
of X^r2^^(X) is equivalent to the inequality 

sTtX^Q^\X) + {1-s)TtY^Q-\Y) (19) 
> Tr {sX + (1 - s)YyQ;l^^,_^)^{sX + {l-s)Y) 

which can be rewritten as 

(20b) 



l-s)Vt-^'J \Y 

s{l-s)\ fx 

sil-s) (i-sr [y 



We now regard A,B,s as fixed and consider the ratio obtained by dividing ()20b|) by 
(j2(Jaj) . Joint convexity holds if this ratio is < 1 for all choices of {X, Y). This is an 

(v Cv) 

optimization problem of the form sup — ^ which is equivalent to the eigenvalue 

{v, Dv) 

problem Cx = XDx when C, D are positive semi-definite. For the operators used 
here, the reduction yields a pair of equations 

^;aV.)b(^^ + (1-^)^) = >^^AiX) (21a) 
n;l_^^,_,)^{sX + {l~s)Y) = Xn-^\Y) (21b) 
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and we want to show that A < 1. Define M = fl^^_^^^__,^^{sX + {l—s)Y) and observe 
that (j21l ) gives three expressions for M since 

M = + (i-.)y) = xn-/ix) = xn-^\Y). (22) 

Applying the appropriate inverse operator to each of these and combining the last 
two, one finds 

sX + {l-s)Y = nsA+(i-s)B{M) (23a) 
X[sX+{l-s)Y] = s^AiM) + {l-s)nB{M) (23b) 

Now Theorem ^ implies 

sTtM^Qa{M) + {1-s)TtM^Qb{M) < Tr Ml"fi,^+(i_,)s(M), (24) 

which then implies 

(1 - X)TtM^[sX + (l-s)y] = (1 - A)Tr Af1'fi,A+(i-s)B(M) < 0. (25) 
It then follows from the fact that QsA+{i-s)b is positive semi-definite that A < 1. 

3.3 From Theorem [2] to Theorem [HI 

Showing that TheoremElimplies the concavity of the map A \—>- Ti e^+^°s^ is elemen- 
tary, but tedious, as it is "merely" a matter of computing derivatives with attention 
to non-commutativity. Let f{x) = Xr e^"'"'"^*^"^"'"^^-' . Then /"(O) < implies that 
rpj,gK+iogA jg concave. One can use power series to verify that ^e^^"^*^ = Q^f{G). 
This can then be combined with p8|) to yield 

/"(O) = TTQA\B)Q^K+io,A[Qj\B)] -2TTe^+^°^^TA{B). (26) 

(Chapter 3 of ,1.3^ gives a clear exposition of the details.) To show that (j2Sl) is 
negative, one can apply the inequality (ji^ to g{x) = G{A + xD, K + xL) with 
G{A,K) = -TrKn^\K). One finds 

g'{0) = -2Tt DTAiK) + 2Tr m^\K) < Tim o\L). (27) 

Now choose K = B,D = e^+'°s^ and L = Vl^k+io^a [^^^(5)] . Then Tr Ln^\L) = 
Ti LQ^^{B). Making the appropriate substitutions in and and ^I7\i yields 

/"(O) < 2 (Tt m^\L) -Ti m^\B)^ = 0. (28) 

Although we have sketched the path from concavity of Tr K'^A^KA^^p to concavity 
of most of the steps are easily seen to be reversible. With a bit more 

effort [7], one can show that the reverse implication also holds. 
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4 Theorem El a direct proof and implications 



4.1 Proof 

Lieb and Ruskai ^U] proved the joint operator convexity of the map {X, A) i— >■ 
XM^^X. Unfortunately, this does not seem to directly imply the joint convexity of 
Tr X'^Q^^{X). However, the following observation allows one to adapt the argument 
in jTU] to give another proof of the joint convexity of Tr X^n^^(X). Let La and Ra 
denote left and right multiplication by A so that La{X) = AX and Ra{X) = XA. 
By expanding in a basis in which A is diagonal, one can verify that 

POD 1 

TrK^n^\K)=Tr / {La + tRA)-\K)-—dt, (29) 
Jo i- + t 



and 



TrATA{K)=Tr / K^La + tRA)-\K)——-dt. 

Jo (1 + ^) 



(30) 



The operator {La + IRa)^^ can be regarded as another non- commutative version 
of multiplication by A'^ . It is positive semi-definite with respect to the Hilbert- 
Schmidt inner product since both La and Ra are self-adjoint and positive semi- 
definite. For example, Tr [La{W)]^X = Tr {AW)^X = TrW^AX = TrW^LAiX) 
and Tr X^La{X) > 0. 

Lemma 4 For each fixedt > the map {A, K) i— >• TrK\LA + tRA)~^{K) is jointly 
convex. 

Proof: Let M,- = {La^ + tRA^)-^/\K^) - {La^ + tRA,)^'\k). Then 
< ^Tr M]Mj 

3 

= J2TrK]{LA,+tRA,r\K,)-Tr{j:^K])A (31) 

3 

- Tr At(^^. K^) + Tr At^ .(L^^. + tL^JA. 

Next, note that {LA,+tRA,){W) = {A,W+tWA,) = {j:jAj)W+tW{J2.A^) 
= L^_4^(VF) + Aj(W^) for any matrix W. Therefore, inserting the choice 

A = {L^^A, + tRY.J,r\j:, K,) in m yields 

Tr (E.i^.yfe.A, +t%,A,)"'(E,/^.) < E,^^K]{LA,+tRA,r\K,). (32) 
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for any t > 0. The joint convexity then follows from the replacement A XjA, 
K ^ XjK with Xj > and Xj = 1. Alternatively, one can simply observe that 
Tr XK\LxA + tRxA)~\XK) = XTi K^{LA + tRA)~\K) and use the first observation 
in Appendix IHl to reduce the joint convexity to (jSH)- QED 

The joint convexity of both TrX^fi^^(X) and Tr AT'^^{K) follow immediately 
from Lemma m and the integral representations above. In particular, inserting (jH^ 
in 1^ yields 

Tr (E,K.)^^^£^^.(E,i^.) < E,TrK]fi:,;(K,) (33) 

which proves Theorem |21 when A = B. The general case then follows as in (fTUI) . 
The choice t = in Lemma E] yields the operator convexity of proved in 

cni. 



4.2 From Theorem [2] to entropy inequalites 

The von Neumann entropy is defined as S{P) = — TrPlogP and the relative 
entropy as H{P,Q) = TrP(logP — logQ) with P,Q positive semi-definite and 
ker((5) C ker(F). (One usually assumes that TrP = TrQ = 1, but this is not 
strictly necessary.) Several special cases of the joint convexity of H{P, Q) follow 
immediately from Theorem |21 

Using (UHl), one can show that ^S(A + xK)\^_^ = — Tr K \og{A), when Tr = 
and that the second derivative satisfies 

-A{A,K) ^ ^S{A + xK)l^^ (34) 

= -Tr Kn^\K) + Tt AT a{K) (35) 

f°° t 
= -Tr / K\LA + tRA)-\K)——-dt (36) 

Jo (1 + ty 

where we used the integral representations (j^^ and ()30|) together with the simple 
fact — (Yip^ = (T+t)2- Lemma m implies that A{A,K) is jointly convex. One 
also has 

H{P + aA,P + bB) =TTA^np\B). (37) 

a=b=0 



dadb 



Now consider density matrices defined on a tensor product of two or three spaces. 
First, suppose that 712 = ^] is block diagonal with A, B in H2 so that 



5(712) = S{A) + SiB) and 72'= Tr 1712 = A + B. Then 

5(712) - 5(72) = 5(A) + SiB) - S{A + B) {3t 
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+ xK 

Now consider the perturbed density matrix 712 = ( „ d 1 r 1 



B + xL 

g{x) = S{A + xK) + S{B + xL) -S[A + B + x{K + L)]. Then 

/(O) = -A(A, K) - A{B, L) + A(A + fi, K + L) < (39) 

because A (A, K) is jointly convex in (A, i^'). This proves the concavity of the map 
712 I— »■ '5(712) — 'S'(72) in the special case of block diagonal matrices. 

The situation above can be compared to P123 = {^^^^ ^t)p" ^ ' ^^^^ case, 

the strong subadditivity inequality in the form S{pi2z) — S{p2z) < S{pi2) — S{p2) is 
precisely the concavity of pi2 ^— > S{pi2) — S{p2) with pi2 = tp[2 + {i — i)Pi2 arbitrary. 
In the previous paragraph, p[2 and p'/j were also required to be block diagonal. 

The joint convexity of the less familiar symmetrized relative entropy follows 
immediately from the joint convexity of TrX^fi^^(X). Using ()29|) one easily finds 

H{P, Q) + H{Q, P) = Tr [{P - " Q)j^] (^0) 

It then follows immediately from Theorem|21with K = P—Q that H{P, Q)+H{Q, P) 
is jointly convex in P, Q. This implies that at least one of H{P, Q) and H{Q, P) is 
jointly convex in P, Q. 

Although Theorem ^ gives a short route to special cases of the joint convexity of 
relative entropy, proving the general case seems to require one of the paths through 
Theorem Q or Theorem 121 Bauman [2j also recognized that there is a connection 
between the concavity of 5'(pi2) — 5'(p2) for block diagonal pi2 and joint convexity 
of Tt X'^Q^^{X) for self-adjoint X; and proved the latter for 2x2 matrices. 



4.3 Further remarks 

The argument used to prove Lemma |3] is remarkably simple^ (and can be applied 
to + tRsj to give a direct proof of the general case of Theorem |21) However, 
the insight needed to rewrite TtK^Q^^[K) in the form (j29j) comes from the very 
powerful relative modular operator formalism originally developed by Araki jT] to 
extend relative entropy to infinite dimensional operator algebras which might not 
even have a trace. The utility of this formalism in finite dimensional situations was 
not realized until much later. For an exposition see Ohya and Petz Lemma 0] 
can also be used to prove the joint convexity of any operator which has an integral 

^In quant-ph/0604206 the argument presented in Section 4.1 is used to give a direct proof of 
the joint convexity of relative entropy. 
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representation as in (j^^ or ()30|) with a positive weight, and the proof presented 
here is similar to one presented in to prove monotonicity of a generahzed form of 
relative entropy and related inequalities within the theory of monotone Riemannian 
metrics developed by Petz |14|. 

In view of connections to the Bures metric and quantum channels, we make a 
few further remarks. Monotone means decreasing under completely positive trace 
preserving maps, and (|29|) is the Riemannian metric associated with the usual rel- 
ative entropy H{P,Q). A large class of such metrics, generalizations of relative 
entropy, and of the operator fi^^ can be obtained by replaced the weight in 
by one coming from a subclass of monotone operator functions. The bound, 

R^' + > Q^' > {Ra + La)-\ (41) 

which holds as an operator inequality is also satisfied by these generalizations of 
n^^, and the upper and lower bounds are special cases of these generalizations. 
This is of interest because any of these Riemannian metrics can be used to define 
a geodesic distance -D(P, Q) between two density matrices. The Bures metric given 

by [-DBures(-P, Q)]^ = 2 [l — Tr (v^QvT^)^''^] is precisely the geodesic distance as- 
sociated with TiX'^{Ra + La)~^{X). The importance of the Bures metric, which 
is closely related to the fideity, in quantum information has been emphasized by 
Uhlmann JH]- Unfortunately, no closed form expression for the geodesic distance 
associated with (piJj) is known. See [21 IH] for bounds and a discussion of the geodesic 
distances for the Wigner-Yanase-Dyson entropies, which are closely related to the 
function considered in Theorem ^ 

Acknowledgment: It is a pleasure to thank Professor M. d'Ariano for the oppor- 
tunity to present lectures during the Quantum Information Processing workshop in 
Pavia, Italy. 

A Maximum modulus principle 

The maximum modulus principle given in most elementary texts, states that the 
modulus of an analytic function can not have a local maximum on a bounded open 
set unless it is a constant. It then follows that a function f{z) which is analytic on 
a bounded set and continuous on its closure achieves the supremum of \f{z) \ on the 
boundary. 

Now consider a function f{z) which is analytic on the strip {2; : < Re^; < 1} 
and continuous and uniformly bounded on its closure {2;:0<Re2;<l}. Let 
M = sup{|/(2;)| : < Re2; < 1}. The maximum modulus principle for a strip 
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[5] says that M is equal to the supremum of \f{z)\ on the boundary, i.e., M = 
sup{|/(2;)| : z = + iy or z = 1 + iy}. 

If the supremum is actually attained at some point z = a + iy, i.e., \f(z)\ = M, 
the result follows easily from the theorem for bounded regions. The possibility of 
a 7^ 0, 1 can be excluded by considering f{z) on the rectangular region Qb = {z : 
< KeZ < 1 and — b < Imz < b} with b > y. Since / cannot have a relative 
maximum on Qf,, one has a contradiction unless / constant. If a = or a = 1, 
then z lies on the boundary of the strip and the assertion holds. To prove the result 
if the supremum is not attained, consider the functions fs{z) = z^'^f{z). Since 
lim \fe{x + iy)\ = 0, the maximum modulus = sup{|/£(2;)| : < Re^ < 1} is 

attained for some z,, = iy or Zf, = 1 + iy. But < and M = limM^ = 

sup{|/(2;)| : z = + iy or z = 1 + iy}. 



B Homogenous concave functions 

Most of the functions considered here are homogenous of degree one, i.e., g{iJ.A) = 
fig{A). In this case concavity is equivalent to superadditivity. To see this observe 
that homogeneity and concavity at | imply 

y{A + B) = g{l[A + B]) > \g{A) + lg{B). (42) 

Conversely, if g is homogenous and superadditive 

g [xA + (1 - x)B] > g{xA) + g[{l - x)B] = xg{A) + (1 - x)g{B) (43) 

One extremely useful property of homogenous concave functions is the following 
derivative inequality [T3] . 

,.^p(X + .B)-,(-4) 

which follows easily from g{A + xB) < g{A) + xg{B). 



C Erratum to |quant-ph/0404126 by M.B. Ruskai 



Since this note appeared, people who have presented the proof in class have re- 
ported some minor errors and omissions in the argument presented in the note. 
The problems and necessary modifications are described in detail below. 
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First, the cyclicity of the trace and procedure described above (0) yield 

\fkiz))\ < \\Ak\\^-''\\C-'\\^-^TT\MC-^''+'y^/^Al+'^C-^''+'y^/^M^\ . (45) 

For y ^ one can not remove the | | because the quantity inside is not positive 
semi-definite. From a pedagogical point of view, the simplest argument may be to 
let G{z) = C-'/^AlC-'/^ and write fkiz) = Tr M^Gkil- z)MGkiz). Then use 
to obtain 

\fk{z)\ < {TiM^Gkil - z)Gkil~z)MY^\TiM^Gkiz)Gk{z)My^\ (46) 

Then, the cyclicity of the trace and proceedure described above (5) yields 

T:tM^Gkiz)Gkiz)M = TTGkiz)Gkiz)MM^ 

< IIC-^^^ pfcf^TrMfM. (47) 

Combining this with a similar estimate for Ti M'^Gk{l — z)Gk{l — 'z)M and taking 
the square root gives the desired result. 

Another alternative is to realize that one only needs a bound independent of z, 
not one with the precise form (jSj). The polar decomposition theorem implies that 
one can write any operator as X = V\X\ with V unitary and |X| = VjOX so that 
|X| = y^X. Using this in (gSD yields 

\fk{z))\ < \\Ak\\\\G-'\\TT\M^V^M\ < ||M||Tr|M|. (48) 

Although the unitary V may depend on z, the last bound on the right is independent 
of z. One could also use ^ to show Tt\M^VW\ < TiM^M = \\M\\l; however, 
the weaker bound above will suffice. 

Second, to bound (jH)) one needs the inequality 

|TrX^r| < (TrX"fX)i/2(TrF"fy)^/2 < 1{Tt X + TiY^Y) (49) 

and ((Tj) should be replaced by 

\fkiO + iy)\ < |(TrM^C*^/2(-7-i/2^^(^-i/2^-iy/2^ ^gg^ 

+ Tr M^G-'y^^G^^/^AkG-'^/^G'y/^M) k = l,2. 

Then (jH)) becomes 

\f{0 + ty)\ < X,\,h{0 + ty)\ + X2\f2i0 + ty)\ 

< \ (Tr M^C*^/2(<7-l/2 ^^^^^ ^ ^^j^^ Q-\l2(j-iyl2^ 

+ llYM'^G"yl^G'^l^{\xAx + A2y42)C"^/^C+^^/^M) 
= \{^xM'^G'yl^G-^yl^M ^'YxM^G-'yl'^G^'yl'^M) = Tr M^M. 

Acknowledgment: The author is grateful to M. d'Ariano and P. Hayden for cor- 
respondence about these issues. 
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